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We argue that calculations in QED at finite temperature are more 
conveniently carried out in the Coulomb gauge, in which only the 
physical photon degrees of freedom play a rol and are thermalized. We 
derive the photon propagator in this gauge for real-time finite tem- 
perature calculations and show that the four-fermion static Coulomb 
interaction that appears in the Lagrangian can be accounted for by 
suitably modifying the photon propagator. The Feynman rules of the 
theory are written in a manifestly covariant form, although they de- 
pend on the velocity 4- vector of the background medium. As a first 
step in showing the consistency and usefulness of this approach, we 
consider the one-loop calculation of the electron self-energy S. It is 
explicitly shown that the divergences that arise from the vacuum con- 
tribution to £ are independent of u^, which implies that the counter 
terms that must be included in the Lagrangian are the same as those 
in the vacuum. 

1 Introduction 

The real-time formulation of Finite Temperature Field Theory (FTFT) is 
very appealing because it can be carried out in a covariant way[Q. This, of 
course, does not come for free. The existence of a preferred frame when the 
background is a material medium, rather than the vacuum, is unavoidable 
and manisfests itself in the dependence of the theory on the vector w M rep- 
resenting the velocity four-vector of the background, which has components 
(1,0) in its own rest frame. Nevertheless, the covariance of the theory thus 
formulated can be exploited in practical calculations, and also allows us to 
deduce very useful results in a general way. For example, the fact that the 
self-energy of a fermion propagating through a medium depends on in 
addition to its momentum k^, implies that the pole of the propagator is no 
longer given by the equation k 2 = 0, which explains why a chiral fermion will 
in general acquire an effective mass0. 

In a theory with fermions and scalars only, the real-time formulation of 
FTFT is straightforward. However, the situation is more complicated in a 
theory with gauge invariance like Quantum Electrodynamics (QED). At fi- 
nite temperature the photon propagator in a covariant gauge (with gauge 
parameter A) can be deduced in several ways]]]. A particularly instruc- 
tive derivation based on a generalization of the Gupta-Bleuler quantization 
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method of the Stueckelberg Lagrangian||, is summarized in Section |2|. As 
discussed there, an essential assumption is implicitly contained in the covari- 
ant formulation; namely, that all the photon degrees of freedom are in thermal 
equilibrium. This is not in accord with the notion that the non-transverse 
photon modes do not appear in the space of states used to define the thermal 
averages. Motivated by a similar reasoning, Landshoff and Rebhan|4| have 
recently proposed a modified expression for the covariant photon propagator, 
in which only the transverse degrees of freedom are in thermal equilibrium. 

We wish to emphasize that the origin of the difficulties lies on the insis- 
tence in formulating the theory in a covariant gauge. In the vacuum, one 
compelling reason for doing so is that, since the theory is Lorentz invari- 
ant, the results of any calculation should ultimately reflect that property. In 
particular, in non-covariant gauges an artificial dependence on some param- 
eter analogous to the vector is introduced, which eventually disappears 
at the end of the calculation of a given physical quantity. Therefore, it re- 
sults convenient to formulate the theory from the beginning in such a way 
that it is manifestly covariant at all stages of a calculation. The situation 
in the medium is different and more complicated. There the presence of the 
background sets up a preferred frame of reference. This is true even if a 
covariant gauge is used, and it can be appropriately accounted for by the 
dependence of physical observables on the vector u^. Therefore, there is no 
apparent advantage in using a covariant gauge in this case. In fact one can 
turn the argument around and argue that, since the theory anyway depends 
on m m , it is natural to use a gauge that also depends on this parameter. 
Morevover, the fact that all the unphysical photon degrees of freedom disap- 
pear in the Coulomb gauge, further suggests that this is a convenient choice 
for implementing QED at finite temperature. 

Motivated by the above reasoning we considered in a previous work[Q 
the calculation and physical interpretation of the the absorptive part of the 
electron self-energy in the Coulomb gauge. Within the formalism of FTFT, 
the damping rate 7 is determined from the imaginary part of the dispersion 
relation of the propagating mode, which is in turn obtained by looking at 
the pole of the corresponding propagator. In Ref. a formula for 7 was 
obtained in terms of the imaginary part of the self-energy, which is valid in 
the physically meaningful situation that 7 and the absorptive part of the 
self-energy are small. It was shown there that the expression for 7 coincides 
with the formula for the total reaction rate V (defined as a combination of 
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probability amplitudes for various processess weighted by appropriate sta- 
tistical factors) ||, provided that: (i) the amplitudes are calculated with the 
properly normalized spinors that satisfy the effective Dirac equation in the 
medium and (ii) the electron self-energy, whose imaginary part determines 
7, is calculated using the Coulomb gauge for the thermal photon propaga- 
tor. On the other hand, if the electron self-energy is calculated with the 
photon progator expressed in a general covariant gauge, then the formulas 
7 and T do not coincide. Even worse, several calculations of 7 carried out 
during the last few years have produced contradictory results, entangled by 
questions of gauge invariance and other problems^. In that sense, the result 
of Ref. || is stimulating and suggests a way for taking a fresh look at the 
subject, guided by what should be a physically reasonable requirement (i.e., 
that 7 = T) which is verified in theories where the ambiguities associated 
with gauge invariance are absent H. 

Eventhough we do not know whether this approach will always yield con- 
sistent results, encouraged by the above considerations we propose to adopt 
the Coulomb gauge as a convenient framework to do calculations in finite 
temperature QED. As a first consistency check, in this work we demostrate 
that the ultraviolet divergences of the electron self-energy do not depend on 
Ufj, and, as a consecuence, can be substracted by means of the same coun- 
terterms as in the vacuum. Therefore, the full Lagrangian, including the 
counterterms, remains Lorentz invariant, despite the fact that we work in a 
non- covariant gauge. 

In Section 0, the formula for the thermal photon propagator in a covariant 
gauge is deduced by adapting the Gupta-Bleuler formalism to the situation 
with a background. This procedure shows in a transparent form the un- 
derlying assumption that the unphysical photon degrees of freedom are in 
thermal equilibrium. In Section |3] we derive the Feynman rules for real-time 
finite temperature calculations in the Coulomb gauge, which were used in 
Ref. M. It is shown that the static Coulomb interaction, which appears at 
the level of the Lagrangian in this gauge, is accounted for by a suitable modi- 
fication of the photon propagator. The photon propagator and the Feynman 
rules are written in a manifestly covariant form, although they depend on 
the velocity 4-vector of the background medium. In Section [| we use 
those Feynman rules to calculate S and demonstrate that the ultraviolet 
divergences are independent of it„- 
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2 Photon propagator in a covariant gauge 



Here we derive the finite-temperature propagator for a massive vector field 
by extending the Gupta-Bleuler quantization method of the Stueckelberg 
lagrangian ||. As explained in the Introduction, this approach, which we 
have not seen in the literature, brings out some of the physical content and 
assumptions that lie behind the covariant formulas for the photon propagator. 
Our starting point is the Stueckelberg Lagrangian, 

L = + l -m 2 A,A» - ^X(d ■ A) 2 . (2.1) 

For A 7^ this Lagrangian admits the correct zero- mass limit, in which we 
are interested. The corresponding equation of motion is 

(d 2 + m 2 )A»-(l-\)d»d-A = 0, (2.2) 

and taking its divergence we obtain 

\^ + Hpj d . A = . (2.3) 

Therefore, for non- vanishing A, d ■ A is a scalar field that satisfies the Klein- 
Gordon Equation with square mass M 2 = m 2 /X. As a consecuence of 
Eq. fl2~3[ ), the field A^ is split into two parts 

A, = Al ~ —^(d ■ A) , (2.4) 

where A^ is a divergenceless spin-1 vector field that satisfies the Klein- 
Gordon equation with mass m. Taking into account the fact that the masses 
of the spin-0 and spin-1 components of A^ are different, its plane wave ex- 
pansion is 



+ 



d 3 k 
2 

d 3 k k 



r d 3 k 3 



i-' 



(2ir) 3 2uJk Tn 



a ko e~ lkx + h.c. , (2.5) 
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where io k and u k are given by 



u k = V k 2 + m 2 



ZJ k = \Jk 2 + M 2 , (2.6) 

and the three space-like orthonormal vectors e At (/c, A) are simultaneously or- 
thogonal to k^ and statisfy 

e(k, A) • e(k, A') = - <5 a ,a' , 

Eau)^(M) = -L*-^). (2.7) 

A=l \ 171 J 

Within this approach the theory is quantized by imposing the indefinite met- 
ric commutation rules 

[a kX ,a* k , xl ] = (2n) 3 2co k S Xt yS^\k-k'), 
Ko,4'o] = - (2tt) 3 2^ (3) (£ - k') , (2.8) 

with all the other commutators being equal to zero. 

The elements of the photon propagator matrix are determined by substi- 
tuting the plane wave expansion of Eq. ( |2.5|) into the following set of relations 

iA lllMU (x-y) = (T(A^x)A„(y))), 

iA 22flu (x - y) = (T (A li (x)A„(y))) , 

iA 12 ^(x-y) = (A u {y)A^(x)) , 

iA 21 ^{x - y) = (A^x)A u (y)), (2.9) 

where the angle brackets denote the thermal average over the states of the 
system and the symbols T and T stand for the time-ordered and anti-time- 
ordered products. The statistical averages of the products of creation and 
annihilations operators are given by 

(a kX a* yx ,) = (2ir) 3 2Lu k 5 x , x ,5 {3 \k-k')(n k + l), 

(a* kX a k , x ,) = (2n) 3 2uj k 5 x , x ,5 {3) (k-k')n k i 

(a k oa* k , ) = ~ (2n) 3 2ZJ k 5 {3 \k - k')(n k + 1) , 

(a* ko a k , ) = - (2n) 3 2uJ k 5^(k - k')n k , (2.10) 
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where 



1 

1 



(2.11) 



with f3 denoting the inverse temperature. 

Using the above the expressions, a straightforward calculation gives 



A^n(A;) 



k^k u 



1 



where 



O 



k 2 — m 2 + ie m 

-2ttzCV [r? 7 (/c) + 9(k ■ u)] , 
-2ni0^ u [r] 7 (k) + 9(-k ■ u)] 

-A;,u(fc), 



knk v 



k 2 — m 2 + ie k 2 — M 2 + 



-^(fc - m 2 ) + 5(r - m 2 ) - - M ) 



(2.12) 



(2.13) 



6* is the step function, and r/ 7 is defined by 



9(k ■ u)tib{x) + 9(— k ■ u)n B (—x) 
1 



Here 



g/3|fc-tt| ^ 

n B (x) = 



1 



e x - 1 

is the boson distribution function written in terms of the variable 

x = f3k ■ u . 



(2.14) 
(2.15) 

(2.16) 



The covariant expression for the photon propagator is obtained by taking 
the zero- mass limit in Eq. ( 2.1 2| ) . Using the relation 



lim — 

m 2 —>0 m 



5{k 2 - m 2 ) - 5{k 2 - M 2 ) 



1 - - 



1\ d 



XJ dk 2 



5{k 2 ) , (2.17) 
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we get 



(*) 

(*) 
A MI/ i 2 (A;) 
(k) 



A-ia/ik) 



k 2 



ie 



2mr] 1 {k)A^ v {k)8(k 2 



A 



(j,v22 



-2ixiA^{k)5{k 2 ) [r, 7 (k)+B{k-u)) , 
-2-KiA^ u (k)5(k 2 ) [j^(jfe) + 0(-Jfe • u)] 
-A^u(fc), 



where 



A 



in' 



9iiv + 1 



(2.18) 



(2.19) 



The above formulas for the photon propagator coincide with the ones 
obtained by the time-path methodp, provided that the derivative of the 
delta function is interpreted according to the prescription given in Eq. ([2.171) . 
Moreover, the derivation presented here (see in particular Eq. fl2.10| )) exhibit 
in a transparent manner the fact that those formulas rest on the unjustified 
assumption that even the unphysical photon degrees of freedom are thermal- 
ized. 

Finally, we notice that the thermal Proca propagator can be obtained 
directly from Eq. Q2.12 ) in the limit A — > 0, with m 2 ^ . This yields 



A^ u (k) = A^(fc) 

A^2i(k) = 
A^ 12 (k) = 



-2ni5(k 2 -m 2 )r](k)^ 



Ativ22{k) 



k 2 — m 2 + ie 
2mA^{k)5{k 2 - m 2 ) [rj(k) + 6(k ■ u)} , 
27ciA^{k)5(k 2 - m 2 ) [rj(k) + 9{-k ■ u)\ 



where 



A 



-9i* + 



kf_i k v 
m 2 



(2.20) 



(2.21) 



3 Photon propagator in the Coulomb gauge 

The plane wave expansion of the free photon field in the Coulomb gauge is 



d 3 k 



(2vr) 3 2^ fc 



£ [a kX e,(k,X)e~ lk - x + h.c.} 



(3.1) 



A=l,2 
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where the polarization vectors e^{k, A) are given by 

e*(M) = MM)), 

with 

e(k, A) • k = , 

and k° = u k = k . We have denoted the vector field by A* r (a;) to indicate 
explicitly that, in this gauge, only the (physical) transverse degrees of free- 
dom are present. The annihilation and creation operators a k \ and a* kX satisfy 
the usual commutation rules appropriate for bosons, 

[a kX ,a* kf y] = (2n) 3 2uj k 5 x ^5^(k-k') . (3.2) 

The photon propagator matrix is determined from the following statistical 
averages 

iAf lflv (x-y) = (T(A t ;(x)A t ;(y))) J 
iA^ u (x-y) = (T(A t ;(x)A t ;(y))), 
iA^ u (x-y) = (AZ{y)A$(x)), 

iA% v {x-y) = (A»{x)A»{y)), (3.3) 

where the symbols T and T have the same meaning as in Eq. ( [2.9|) . The 
calculation of the propagator is identical to the one presented in the previous 
section, with some obvious modifications. Thus, substituting the plane wave 
expansion of A^ in Eq. (|3.3| ) and following steps similar to those that lead 
to Eq. ( j2.12j ), in the present case we obtain 

= (-R,u)A ab (k) , (3.4) 



where 



A u (k) = ■^--2mr ]l 5(k 2 ), 

A 22 {k) = jf^— - 2ixir)^{k 2 ) , 

A ia (ife) = -27Ti5(k 2 )[7]^ + e(-k-u)}, 

A 21 (k) = -2^5(^)1^ + 6(k-u)}, (3.5) 
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with r/y defined in Eq. ( |2.14| ). The tensor R^ v is given by 



(3.6) 



A=l,2 



and its explicit expression in terms of and fc M is 



R, 




V 



00 



J , , , . h 
^{u^k v + k^Uy) + —Uf/u 



(3.7) 



where 



k ■ u , 



Vto 2 — k 2 



are the energy and the magnitude of the 3-momentum in the rest frame of 
the medium. It should be noted that the term in R^ v depending on u^Uy 
disappears from the background-dependent part of the propagator because it 
is proportional to k 2 5(k 2 ). That term also disappears from the background- 
independent part for the reason that we explain below. 
In the Coulomb gauge, the interaction Lagrangian is 



with p = jo- Therefore, in addition to the contribution involving the pho- 
ton propagator between electron lines, there appears a static four-fermion 
interaction that in principle must be taken into account in any calculation. 
However, it is possible to absorb the effect of the four-fermion interaction into 
the photon propagator by the following device. The interaction in Eq. (|3.8| ) 
can be thought of as being produced by the exchange of a scalar particle 4> 
with the following propagator 




(3.8) 



where e is the electron charge and, as usual, 



(3.9) 



iA^(k) 



(3.10) 



K, 2 — ie 
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and an interaction Lagrangian with the electron of the form 

= -ep<f) (3.11) 

Therefore, instead of the interaction Lagrangian given in Eq. ( |3.8| ), for cal- 
culations we can use 

Lim = ~ep<P + eA tr ■ j , (3.12) 

where, in order to reproduce the effect of the four-fermion interaction, must 
be assigned the following propagator 



k 2 — it 



K + l€ 

A { $ = A ( i> = 0. (3.13) 



Eq. ( 3.12 ) can be written in a compact form by introducing the field 

A,=A^ + u^. (3.14) 
Then, instead of Eq. ( |3.12| ) and the two separate propagators A^ b and 



A (0)afe> we use 

L int = -ej^A^ (3.15) 
together with the combined propagator 

A abflv = + u,u u A[f . (3.16) 
From the formulas in Eq. ( |3.4j ) and Eq. (|3.13| ) we finally obtain 

A abfll/ {k) = (S^)A ab {k) , (3.17) 

where S^ u is given by 

Snu = 9uv H ^knk v ^{Unky + knUy) (3.18) 

It is useful to observe that, at to = k, = and therefore, according 
to Eq. ( p.6|) , it follows that 



S ^L=K = ~ Yl e M (MMM) 



A=l,2 



(3.19) 



To summarize, in the Coulomb gauge the photon propagator can be taken 
to be A ab ^ v as given in Eq. (|3.17| ), where the elements A ab are specified in 
Eq. (|3~5|) , and with the interaction Lagrangian of Eq. Q3.15| ). 
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4 Divergences of E 



Now, we turn our attention to the subject of the ultraviolet divergences of 
the real part of the electron self-energy and the way to dispose of them in 
calculations done within the Coulomb gauge. Since the background depen- 
dent parts are finite, we only have to worry about the vacuum term, which 
is given by 

-iE<°> = He) 2 / ^riSP(p + k)^iA^(k) (4.1) 



where 



and 



SP = *±* (4.2) 



(p + k) 



k^ky k ■ u . i.i \ 



with k, 2 = (k.u) 2 — k 2 . The contribution produced by the term proportional 
to g^ v is the usual one calculated in the Feyman gauge, and it does not 
depend on u^. Our aim is to show that the contributions from the other two 
terms of the photon propagator can be split into a finite part that depends 
on Ufj,, and a divergent part that is independent of it. To do that we consider 
the quantity 

1 = J ^"tVCp + fc)aM» - k ■ «M* + K^)] > ( 4 - 4 ) 

where 

D = k 2 (p + k) 2 K 2 . (4.5) 

Using the identity 

VYY = \(9 m 9 uX - 9" V 9 aX + s'V" + ^V)7A , (4.6) 



the integral in Eq. ( |4.4| ) can be expressed as 

1 = \j^ X {H(P + k) 2 -p 2 -2{k-u)(p-u)]+p x [2{k-u) 2 -k 2 } 

- 2u x (k ■ u)[(p + kf - p 2 - p ■ k]} (4.7) 
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The terms with the factor (p + k) 2 produce integrals of the form 



1 



<?k 



2^2 ' 



k 2 K 



which vanish upon symmetric integration, while those proportional to p 2 are 
finite as is easily seen by the ordinary power counting argument. The rest 
can be rewritten in the following way 

/ = \l^(p, u) {m^u u - + 2u^[ Pu i - 7 „(p • u)]} (4.8) 

with 

I^p,u)= jd'k^. (4.9) 

By considering the Taylor expansion of the integrand with respect to p, we 
see that the difference I^ip, u) —1^(0, u) is a finite quantity. Thus, in order 
to examine the divergent contribution to /, in Eq. ( |4.8| ) I^ v (p,u) can be 
replaced by 7 Ml/ (0, u), which has the structure 

I M „(0, u) = ag^ u + bu^Uu . (4.10) 

The coefficients a and h are scalar quantities and can only depend on u 2 , 
which is equal to one. Therefore, although they are infinite, a and b do not 
depend on u^. From Eq. ( |4.10| ), we have 

1^(0, oc u u , (4.11) 

which implies that, when contracted with 1^(0, u)u^ the term inside the 
second square bracket in Eq. (|4.8|) vanishes identically. Then, we are left 
only with the divergences arising from the terms in the first square bracket 
of this equation. They do not vanish in general, but are proportional to 

both of which are independent of u M as a consequence of Eq. ( [4.10 ). This 
complete the proof of the statement we made in the Introduction that, al 
least at the one-loop level, the divergences of the electron-self energy are 
independet of the velocity four-vector of the medium. 
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5 Conclussions 



The use of a non-covariant gauge, such as the Coulomb gauge, is inconvenient 
in the vacuum because it requires that the quantization be carried out in a 
particular frame. In contrast, the presence of the medium defines a preferred 
frame which can be exploited to use a gauge in which the unphysical degrees 
of freedom disappear alongwith the associated question of whether they have 
a thermal distribution or not. 

In this article we examined in detail the possible dependence on of 
the divergent contributions to the electron self-energy calculated within the 
Coulomb gauge formulation of QED at finite temperature. We explicitly 
showed that no such dependence ocurrs at the one-loop level and, therefore, 
the renormalized Lagrangian of the theory preserves its Lorentz invariant 
structure, despite the fact that we are using a non-covariant gauge. Taking 
into account the difficulties encountered in the calculations of the fermion 
damping in the covariant gauges, the technical result of this paper reinforces 
the conclussions of Ref. || emphasizing the use of the Coulomb gauge as an 
appropriate choice to carry out calculations in finite temperature QED. 
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